Evidence for long-lived quasiparticles trapped in superconducting point contacts 
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We have observed that the supercurrent across phase-biased, highly transmitting atomic size 
contacts is strongly reduced within a broad phase interval around n. We attribute this effect to 
quasiparticle trapping in one of the discrete sub-gap Andreev bound states formed at the contact. 
Trapping occurs essentially when the Andreev energy is smaller than half the superconducting gap 
A, a situation in which the lifetime of trapped quasiparticles is found to exceed 100 [is. The origin 
of this sharp energy threshold is presently not understood. 

PACS numbers: 74.45. +c,74.50.+r,73.23.-b 



Both theory and experiment indicate that the number 
of quasiparticles in superconductors decreases exponen- 
tially as the temperature is lowered, while their recombi- 
nation time increases [1, 2[. This slow dynamics is an im- 
portant ingredient in non-equilibrium superconductivity 
and allows for the design of high-performance devices like 
single photon detectors for astrophysical applications [3] . 
However, recent developments on microwave resonators 
[4] and Josephson qubits [5] show that at very low tem- 
peratures residual non-equilibrium quasiparticles set a 
limit to the proper functioning of these devices. More 
drastically, a single quasiparticle can determine the re- 
sponse of single-Cooper pair devices [6] containing small 
superconducting islands in which the parity of the total 
number of electrons actually matters. The trapping of a 
single quasiparticle in such a superconducting island has 
been dubbed "poisoning" [7], as it inhibits the behavior 
expected in the ground state of the system. Remark- 
ably, it has been argued [8] that quasiparticle trapping 
could also occur in the discrete Andreev bound states 
[9] formed at sub-gap energies in a constriction between 
two superconductors, a system containing no island at 
all. We demonstrate this phenomenon with an experi- 
ment on atomic size constrictions, where the trapping of 
a single quasiparticle is revealed by the full suppression 
of the macroscopic supercurrent through its well trans- 
mitted conduction channels. We also show that, as an- 
ticipated [8], trapped quasiparticles are long-lived, with 
time scales up to hundreds of fis. 

We use micro-fabricated mechanically controllable 
break junctions [10] to obtain aluminum atomic point 
contacts embedded in an on-chip circuit, sketched in 
Fig. 1(c) [11]. The circuit allows measuring for each 
atomic contact both the current-voltage characteristic, 
from which one determines precisely the ensemble {t{\ of 
the transmissions of its conduction channels [12], and its 
current-phase relation [13]. In order to go reversibly from 
voltage to phase biasing, the atomic contact is placed in 
parallel with a Josephson tunnel junction (critical current 
Io ~ 554 n A much larger than the typical critical current 
of a one-atom aluminum contact ~ 50 nA) to form an 
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Figure 1. (Color online) (a) Short one-channel constric- 
tion between two superconducting electrodes (phase differ- 
ence 5 = 5l — Sr). (b) Excitation spectrum: besides the 
usual continuum of states above the energy gap A, that ex- 
tends all across the structure, there is at the constriction a 
discrete Andreev spin-degenerate doublet with energy Ea(S) 
above the ground state, where quasiparticles can get trapped. 
(c) Schematic setup: an atomic point contact (red triangles) 
forms a SQUID with a Josephson junction (green checked 
box). Phases S and 7 across contact and junction are linked 
by the flux <j> threading the loop. The SQUID is biased by 
a voltage source Vb in series with a resistance Rt = 200 Q. 
The current I is measured from the voltage drop across Rb- 
An on-chip antenna is used to apply fast flux pulses to the 
SQUID. It is represented on the right hand side as an induc- 
tor current-biased with a source If. (d) Four possible config- 
urations of a one-channel constriction: ground state (energy 
—Ea), Andreev doublet empty; two odd configurations, with 
zero energy and definite spin ±1/2, one quasiparticle added 
to the contact; last configuration corresponds to spin-singlet 
double excitation (energy +Ea). 



asymmetric dc SQUID. An on-chip antenna allows ap- 
plying fast flux pulses through the SQUID loop and a 
superconducting coil is used to apply a dc flux. 

In the usual semi-conductor representation, there is 
just one pair of Andreev bound states in a short 
one-channel constriction, with energies ±Ea{S,t) = 



± A y 1 — t sin 2 (5 / 2) determined by the channel trans- 
mission t and the phase difference S across it [8, 11, 14, 
15]. In the ground state, only the Andreev bound state at 
negative energy is occupied, leading to a phase dependent 
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term — Ea(5, t) in the total energy, and a supercurrent 
—Ia = —<-Po 1 dEA/dS, with ipo = K/2e. The two Andreev 
bound states give rise to the excitation spectrum shown 
in Fig. 1(b), with a discrete spin-degenerate doublet, lo- 
calized at the constriction, at an energy Ea < A above 
the ground state. The four lowest-lying configurations 
of the system are built from this doublet. Above the 
ground state, there are two "odd" configurations (spin 
1/2) with a single excitation of the doublet at Ea, i.e. 
with a quasiparticle trapped in the constriction. In this 
case the global energy is zero, i.e. phase independent, 
and the total supercurrent is zero [15]. Finally, there 
is another spin-singlet configuration with a double exci- 
tation, which carries a supercurrent +Ia exactly oppo- 
site to the one in the ground configuration. Hence, the 
supercurrent through the constriction is a probe of the 
configuration of the system. 

In our experiment, the supercurrent through the 
atomic contact is accessed through measurements of the 
"switching current" of the SQUID, which is the bias cur- 
rent at which the whole device switches from the super- 
current branch (V = 0) to a finite voltage state. Because 
of the large asymmetry, the SQUID switching current is 
only slightly modulated around that of the junction by 
the applied flux <f>. The modulation corresponds essen- 
tially to the current-phase relation of the atomic contact 
[13]. As the SQUID loop is small, the phase 7 across the 
Josephson junction, the phase 8 across the atomic con- 
tact, and the phase ip = 4>l<-Po related to the external flux 
<f>, are linked [11] through ip = 5 — 7. To measure the 
switching current of the SQUID, current pulses of vari- 
able normalized height s = I /Iq are applied through the 
bias line, while monitoring the voltage across the SQUID. 
To ensure that the measurements are statistically inde- 
pendent, additional short prepulses that force switching 
are applied before each one of them [11] (top left inset of 
Fig. 2). The switching probability P sw (s) is obtained as 
the ratio of the number of measured voltage pulses to the 
number of bias pulses (typically 10 4 ). In Fig. 2, we show 
P sw (s,(p) measured at 30 mK on one particular SQUID 
({rj ={0.994,0.10,0.10}). For most flux values, we ob- 
serve the generic behavior for Josephson junctions and 
SQUIDs, i.e. a sharp variation of the probability from 
to 1 as the pulse height is increased (lower left inset of 
Fig. 2). However, in a broad flux range 0.7ir < ip < l.ln 
around n. the behavior is completely unusual: P sw (s) in- 
creases in two steps and displays an intermediate plateau 
(top right inset of Fig. 2). 

Precise comparison between experiment and theory is 
performed using an extension [11, 13] of the well-known 
model describing the switching of a Josephson junction 
as the thermal escape of a particle over a potential bar- 
rier [16]. For our SQUIDs, the potential is dominated 
by the Josephson energy of the junction but contains a 
small contribution X^ c i^4(7 + V^i) from the atomic 
contact, which depends on {7^} and on the configuration 
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Figure 2. (Color online) Color plot of measured switch- 
ing probability P 3m (s,<p) for SQUID with contact transmis- 
sions {0.994,0.10,0.10}. Top left inset: Measurement proto- 
col. Short prepulses ensure same initial conditions before each 
measurement pulse of height / = slo and duration t p =lus 
(the subsequent lower plateau holds the voltage to facilitate 
detection). Delay between prepulse and measurement pulse 
is here At = 2/xs. Main Panel: Black curves: theoretical pre- 
dictions (solutions of P sm (s,ip) = 0.5) for pristine (solid line) 
and for poisoned contact (dashed line). Insets: measured 
Psm(s) (red dots) at fixed flux. Lower left: p — (green line 
in main panel). Upper right: ip — 0.88-7T (cyan line in main 
panel). In both insets Ps W (s) (Pj tl) (s)) ) the solid (dashed) line 
is the theory for the pristine (poisoned) contact. In the upper 
right inset, the intermediate line (gray) is a fit of the data 
with the linear combination of Eq. (1) with p — 0.36. 



of its Andreev levels (ci = — 1 if channel i is in its ground 
state, Ci = 1 in excited singlet, and Ci = in an odd con- 
figuration). The predictions for P sw (s) are shown as lines 
in the insets of Fig. 2. Whereas the data taken at (p = 
are well fitted by theory with all channels in the ground 
state, those taken at <p = 0.887T are not. However, they 
can be very precisely accounted for by the weighted sum 
of the theoretical curves Pg W (s) and Pg W { s ) correspond- 
ing, respectively, to the "pristine contact" (i.e. with all 
channels in their ground configuration), and to the "poi- 
soned" contact (i.e. with its most transmitted channel 
in an odd configuration). This is the case in the whole 
flux region where the measured switching curves have an 
intermediate plateau, showing that 

P SW (S, <p) = (l- P M)P1(S, ip) + p fo>) PL (S, if) . (1) 

The function 1 — p ((p) describes the height of the in- 
termediate plateau in P sw (s). A similar analysis was 
performed for other SQUIDs formed with atomic con- 
tacts having one highly transmitted channel [11] and in 
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Figure 3. (Color online) Relaxation time T\ (blue points) 
and asymptotic poisoning probability j>oo(red points) as a 
function of the phase 5 imposed for a time At between pre- 
pulse and measurement pulse. Data taken on SQUID with 
atomic contact {0.994,0.10,0.10} at 30 mK. Measurements at 
\S — 7r > 0.5-7T show very fast relaxation that could not be 
resolved reliably in our setup. Left inset: 1 — p is the height 
of the intermediate plateau in P sw (s). Dashed line: P S m{s) 
found for At ;§> lOOus. Right inset: Typical time evolution of 
p, from where Ti and Poo are extracted. 



all cases, P sw (s) shows a plateau delimited by the predic- 
tions for the pristine and the poisoned contact in a broad 
phase range around n. The fact that the data are pre- 
cisely accounted for by this linear combination induces 
us to interpret the coefficient p as the poisoning proba- 
bility, i.e. the probability for the atomic contact to have 
a quasiparticle trapped in its most transmitting channel. 

We have found that at fixed s and ip, the poisoning 
probability p depends exponentially on the delay At be- 
tween the prepulse and the measurement pulse (Fig. 3 
right inset): a fit of the form p(At) — p^ + (po — 
Poo) exp(— At/Ti) gives the initial poisoning just after the 
prepulse po , the asymptotic value at long times p^ , and 
the relaxation time T\ . To obtain a meaningful measure- 
ment of the phase dependence of the relaxation, we had 
to implement a refined protocol [11] involving flux pulses 
applied through the fast flux line within the time interval 
At. It allows probing the relaxation from a fixed po, with 
measurement occuring always at the same flux, the only 
adjustable parameter being the phase S during the wait- 
ing time. Both p^ (<$) and T± (5), measured at 30 mK, are 
shown in Fig. 3 for the same SQUID as in Fig. 2. Their 
phase dependence is symmetric and peaked at S = n, 
where the relaxation is the slowest and p^ the largest. 
A rapid decay of T\ by almost two orders of magnitude 
and a drop of p^o to are observed at \S — ir\ ~ 0.3-7T. The 
overall shape of both p^ (5) and T\ ($) remains very sim- 
ilar when temperature is varied [11]. The relaxation time 
T\ falls rapidly with temperature, and becomes too short 
to be measured above 250 mK. Similar data [11] taken on 
a variety of atomic contacts show that the phase interval 
in which poisoning occurs reduces when the transmis- 
sion of the most transmitting channel diminishes. For 



channels with all transmissions smaller than 0.7, the poi- 
soning probability p was too small to be measured. An 
important observation is that when two switching pre- 
pulses are applied (instead of a single one) with more 
than lus delay between them, the first one has no effect, 
which indicates that quasiparticles created by switching 
matter only during 1 us. This experimental observation, 
plus the fact that diffusion is expected to efficiently drain 
away from the constriction the quasiparticles created by 
the prepulse [11], allow us to conclude that the resid- 
ual quasiparticle density is constant during the poison- 
ing probability relaxation, and that it does not originate 
from the switching pulses (contrary to po). The fact that 
Poo 7^ proves that quasiparticles are present in the con- 
tinuum in the steady state, as found in other experiments 
[4, 17, 18]. 

When quasiparticles jump between the Andreev states 
and the states in the continuum, transitions arise among 
the four configurations accessible to the Andreev dou- 
blet, as shown in the inset of Fig. 4. All the microscopic 
processes involved are in principle rather slow because 
they either require energy absorption or the presence of 
quasiparticles in the continuum [8[. We define the rates 
Tin and r out corresponding to an increase or a decrease of 
the number of quasiparticles in the contact. Because we 
see no trace of the state with a double excitation neither 
in these data, nor in preliminary spectroscopic measure- 
ments, we assume that the relaxation rate to the 
ground state is much larger than Fi n and r out . From a 
simple master equation [11] for the population 1— p of the 
ground state, and p/2 of each of the odd configurations, 
one gets 



T, = 



1 



2r, 



r f! 



3r,, 



3r in 



(2) 



and the flux dependence of Ti n and r out can then be 
extracted from the data. Then, instead of plotting the 
results as a function of the applied flux 6 like in Fig. 3, 
we choose for the x-axis the Andreev energy Ea (5) of 
the most transmitting channel. The results for five con- 
tacts are shown in Fig. 4, together with the relaxation 
time T\ and asymptotic poisoning probability p^. Al- 
though the Andreev energy is clearly not the only rele- 
vant parameter, the rates for all contacts roughly coin- 
cide. The most apparent differences are in the asymp- 
totic poisoning probability p^ which, for a given Ea, di- 
minishes when the transmission of the most transmitting 
channel increases. Two distinct regimes are evidenced 
in Fig. 4: when Ea/A > 0.5, the relaxation time is 
very short and the asymptotic poisoning is negligible. 
In terms of rates, F; n is smaller than r out by 2 to 3 or- 
ders of magnitude. In contrast, when the Andreev en- 
ergy lies deep in the gap (Ea/A < 0.5), relaxation is 
much slower and the asymptotic poisoning probability 
becomes sizable. This regime corresponds to a smaller 
ratio r out /ri n . The sharp threshold at Ea/A ~ 0.5, 
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Figure 4. (Color online) Relaxation data for five different 
atomic contacts (transmissions are given in panel (b)): (a) 
relaxation time Ti, (b) asymptotic poisoning probability p^, 
(c) rates Ti n and r out as a function of normalized Andreev 
state energy Ea/& of most transmitting channel. There is 
a sharp threshold at Sa/A ss 0.5 for all contacts. Minimal 
value of -Ea/A is \Jl — r (0.08 for the black points, 0.5 for the 
orange points) and is reached when 5 — n, as shown in (d). 
Inset of panel (c): Rates Ti n (resp. r out ) are for processes 
increasing (resp. decreasing) the number of quasiparticles in 
contact. The relaxation rate from the excited singlet to the 
ground state is assumed to be much faster than all other rates. 
The occupation of all four configurations is given in italic 
letters: 1 — p for the ground state, p/2 for each of the odd 
configurations, and for the excited state. 



where r out drops by two orders of magnitude, is observed 
for all the measured contacts with highest transmission 
above 0.74. Furthermore, no poisoning was observed in 
contacts in which all channels had transmissions below 
0.74, the Andreev state energy then being always larger 
than AVI - 0.74 ~ 0.5 A. In contrast, we have found 
that in contacts with more than one highly transmitting 
channel, poisoning can affect several channels at once 
[11]. Presently, we do not have an explanation for the en- 
ergy dependence of the rates: the mechanisms commonly 
used to describe quasiparticle dynamics in superconduc- 
tors (recombination, phonon emission and absorption) do 
not lead to such a sharp threshold at energy 0.5 A. 

Let us mention that it is possible to untrap quasi- 
particles. For example, we have implemented an effi- 



cient "antidote" protocol [11] based on dc flux pulses that 
bring the Andreev states at the gap edge from where 
the trapped quasiparticle can diffuse away into the elec- 
trodes. Furthermore, it is possible to avoid altogether 
poisoning: when the large scale on-chip wires connecting 
to the SQUID are made out of either a normal metal [13] 
or a superconductor with a lower gap than the device 
[11], they act as good quasiparticle traps and poisoning 
is never observed. 

To conclude, we have performed the first observation 
and characterization of single quasiparticles trapped in 
Andreev bound states. The long lifetimes that we have 
measured open the way to individual spin manipula- 
tion and to superconducting spin qubits [8]. Moreover 
the complete suppression of the macroscopic supercur- 
rent when a single quasiparticle is trapped shows that 
a superconducting quantum point contact can be seen 
as a very efficient quasiparticle detector. Finally, let us 
mention that quasiparticle trapping, which is likely to 
be a generic phenomenon in superconducting weak links, 
could be detrimental in some situations. It could be the 
case for experiments proposed to detect "Majorana bound 
states" in condensed matter systems [19] since their topo- 
logical protection relies on parity conservation. 
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We first discuss the various representations of the spectrum of Andreev bound states, then the 
sample design and fabrication, and the switching current measurements. The details of the poisoning 
dynamics measurements are then reported, followed by a description of a procedure, nicknamed 
"poisoning antidote", to untrap quasiparticles from the point contact. A discussion of the various 
processes at play follows, with a presentation of the results in terms of poisoning and unpoisoning 
rates as a function of the energy of the Andreev states. Details on the dynamics of our circuit are 
discussed in appendix A. Appendix B focuses on the calculation of the relaxation rate from the 
spin-singlet double excitation. 

PACS numbers: 



I. ANDREEV BOUND STATES 
REPRESENTATIONS 

The widely used picture for Andreev bound states 
(ABS) in a short single channel conductor is the semicon- 
ductor representation 2,3 shown in Fig. 1, assuming spin 
degeneracy (continuum states at energies smaller than 
—A and larger than A are not shown). The lines corre- 
spond to the energies of the available energy levels. In 
the ground state, all levels at negative energies are occu- 
pied. The total energy has a single phase-dependent term 
—Ea(5). arising from the lowest ABS. When adding a 
single quasiparticle to the system (quasihole or quasielec- 
tron), one accesses the two odd configurations, which cor- 
respond to the ABS being both empty or both occupied. 
Both configurations have zero energy (— Ea+Ea or 0+0), 
in absence of magnetic field or spin-orbit coupling 4 ' 5 . Fi- 
nally, adding two quasiparticles or exciting the system 
with photons at energy 2Ea gives access to the spin- 
singlet double-excitation configuration, represented here 
with the lowest ABS being empty and the top one being 
occupied. The excitation spectrum shown in Fig. Id of 
the Letter directly follows: at a given phase 5, two pos- 
sible single particle excitations have an energy Ea above 
the ground state energy; the cost of the double excitation 
is 2E A . 

An alternative representation, used in Ref. 5 shows the 
energy of the various configurations as a function of the 
phase difference. Without spin-dependent interactions, 
one obtains Fig. 2. 

As discussed in the main paper, the vanishing of the 
supercurrent is interpreted as arising from poisoning of 
the Andreev doublet by a single quasiparticle, i.e. a 
transition from the ground state|0) to an odd configu- 
ration 1 1 X) or 1 1 I) . We think that the following de- 
scription, although not rigorous, helps clarifying the na- 
ture of the Andreev states and their spin structure: in 
the ground state of the system, a Cooper pair state 
(spin singlet state) is formed at the contact with an en- 
ergy —Ea (6) £ [— A,0], through Andreev reflections in 
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Figure 1: Phase (8) dependence of the Andreev bound states 
energies ±Ea in a short transport channel of transmission r. 
The energy A is the superconducting gap. 

the superconductors on both sides. Mimicking the BCS 
theory in second quantization formalism 2 , and inspired 
by Ref. [6,7], the superconducting ground state wave- 
function at the contact can be thought of as 

|*o) = |*cont) ® (ua + v A c\ t c A ^) |0) , (1) 

with lucent) describing the states with energy lying be- 
low —A (continuum states) and |0) is the vacuum. We 
have singled-out from the BCS wave-function the state 
"A" corresponding to the lowest energy Andreev state. 
The creation operators c A ^ j_ create an electron at the 
Fermi level in the channel, with either spin. Hence, the 
coefficients ua and va differ only by a phase, and 

|*o) = l* cont > ® ^ (l + e'MA) 10) • (2) 

Thus, the Andreev energy —Ea can be seen as the con- 
densation energy of the "Andreev Cooper pair". 

The addition of one quasiparticle to the system at 
the lowest possible energy is described by one of the 
two Bogoliubov operators 7^ = u A c A ^ — v a c M or 
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Figure 2: Energy of the different configurations accessible 
to a single channel. The energy of the ground configura- 
tion, labeled |0), is that of the lowest Andreev level — Ea 
of Fig. 1. Two odd configurations, labeled |lf) and |1 4), 
with one quasiparticle added to the ground state, have zero 
energy and a definite spin 5 . In the configuration |2), only the 
excited Andreev state +Ea is occupied. 

<Al = u *a c ai + v A c At- which lead to the "odd" states 

|*if)=7jitl*o> = |*co»t>®c] 4t |0) (3) 

and 

l*U>=7jul*o> = 1*^)8^10), (4) 

both having a single electron occupying with certainty 
the state at the Fermi energy. The energy of the odd 
states is higher than the ground state energy by the con- 
densation energy Ea of the Cooper pair in the Andreev 
bound state. For completion, the upper Andreev state 
|^2) is obtained when two quasiparticles are added to 
the ground state, at an energy cost : 

l*2>=7M t l*o> = -7^7^1*0) 

/ t t \ ( 5 ) 

= |*cont) ® yu* A c\^ Ai -V A j |0) . 

Apart from a non-physical global phase, this state can be 
rewritten as |* cont ) ® 75 (~! + e i5A c Af c Ai) 1°) > which 
resembles the ground state \^q): it contains an excited 
Cooper pair (excited singlet). As a matter of fact, the 
present experiment was initially designed to perform the 
spectroscopy of the transition between \*S? ) an d 1*2} , a 
goal that we haven't reached yet. 

II. SAMPLE FABRICATION AND 
MEASUREMENT SETUP 

The samples were fabricated on a polished, 500 jim- 
thick Kapton substrate, which is insulating and elastic. 
A suspended bridge is fabricated by e-beam lithography 
of a thin constriction in an aluminum film, and etching 



of 1 um of an underlying polyimide layer. Optical micro- 
graphs of the sample at various scales, a SEM micrograph 
of the core of the sample and a photograph of the break- 
junction mount used to obtain atomic contacts are shown 
in Fig. 3. The bridge is placed in a superconducting loop 
containing also a 2.8 nm 2 Josephson junction fabricated 
at the same time by shadow evaporation (70 nm Al, ox- 
idation in an Ar — O2 (85%-15%) mixture at 18 mbar 
for 5 min, 54 nm Al), hence forming a SQUID 8 ' 9 . The 
SQUID is connected with thin (124 nm), narrow (0.9 pi) 
and long (0.4 and 0.8 mm) Al wires corresponding to a 
total inductance L in the nH range. At a larger scale, the 
connecting wires overlap over 0.01 mm 2 a floating under- 
lying 30 nm-thick aluminum electrode, which was covered 
with five heavily oxidized 1.5 nm-thick aluminum layers, 
forming a capacitor C ~ 65 pF, hence shunting the con- 
necting lines at high frequency. The resonance frequency 
(27rV LC)~ x ~ 0.5 GHz and the losses of the correspond- 
ing tank circuit (see Fig. 4) , modelled by series resistance 
r ~ 0.5 SI, were measured by microwave reflectometry in 
a separate experiment with the same sample. From the 
capacitor, a 5017 coplanar waveguide makes the connec- 
tion towards a mm-size connection pad placed at the edge 
of the sample. 

The sample was cooled in a dilution refrigerator op- 
erated down to 25 mK. The substrate is clamped at one 
of its edges, between a small copper plate placed under- 
neath and the half-cylindrical central pins of two SMA 
launchers, which connect to the bias and fast flux lines 11 
(see Fig. 3(e)). Connections of the ground plane to the 
sample holder are achieved in the same manner 12 . The 
atomic contacts are obtained by bending the sample with 
a pusher placed on the side opposite to the contact pads. 
The contacts formed with this setup are very stable. For 
example, the contact discussed in the Letter, labeled AC1 
in the following, was kept unchanged during 2 months, 
till we decided to form another one. A shielded supercon- 
ducting coil (which had been withdrawn before taking the 
picture in Fig. 3(e)) is placed a mm above the SQUID, 
for dc flux biasing. Bias lines connected to the SQUID 
and to the antenna are coaxial lines heavily attenuated 
at various stages of the cryostat, with total attenuation 
of 55 and 30dB, respectively, in order to damp the cur- 
rent noise of the components placed at higher temper- 
ature. The current and voltage measurement lines are 
twisted pairs, equipped with micro-fabricated microwave 
filters 13 . 



III. DETERMINATION OF THE 
TRANSMISSION COEFFICIENTS 

The current-voltage characteristic of the Josephson 
junction alone, taken after opening completely the 
atomic contact, is shown with black open symbols in 
Fig. 5. As observed in other experiments with similar 
junctions 9,11 ' 14 , some current is found at sub-gap volt- 
ages, here for |V| < 200 /ieV. The large scale resistance 
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Figure 3: (a) Large scale photograph of the sample, showing the end of the fast flux (FF) and current bias coplanar lines. 
The bias line and the ground plane form planar capacitors with a common, floating and superconducting rectangular electrode 
placed below, implementing the capacitor C. From there the connections to the SQUID are made with 0.9 um-wide Al wires 
which are not visible at this scale, and have been redrawn for clarity. They appear on the photograph (b) as three dark lines 
with regularly spaced bright spots, one connected to the bias line, the two others to the ground plane. The bright spots are 
wider pillars that hold on the substrate since they are large enough not to be completely freed during the etching step. The fast 
flux coplanar line ends with a short close to the SQUID, the current in the upper half of the short creating a magnetic flux in 
the loop. Photograph (c) shows the SQUID loop, with the Josephson junction on the right and the suspended bridge where the 
atomic contacts are formed on the left. Five small square gold electrodes, intended as quasiparticle traps 10 , are barely visible 
through the Aluminum layers, (d) SEM micrograph of the SQUID, seen under an angle: two angle evaporation of aluminum 
define the superconducting loop with a tunnel junction on the right arm, and a suspended micro-bridge on the left arm, which 
is broken at low temperature to form atomic contacts, (e) View of the sample holder, with a bent sample. The sample is 
clamped between the launchers of two SMA connectors, which are visible on the right-hand side, and a small metallic plate 
hold with two screws at the bottom. The sample is bent with a brass blade held by a rod moving vertically (on the left). For 
this photograph, the bending was exagerated, and the superconducting coil placed immediately above the sample was removed. 




Figure 4: Full schematic of the sample. The SQUID is biased through an on-chip LC circuit (L ~ 1.6 nH, C ~ 65 pF). The 
losses in the capacitor are modelled by series resistance r ~ 0.5 Q. The voltage V measured across the full on-chip circuit 
corresponds, at low frequency, to the voltage drop across the SQUID. 



is Rjj — 550 SI. The same figure presents with red solid 
symbols the I — V characteristic of a SQUID with an 
atomic contact: as compared to the previous curve, the 
large scale conductance is slightly increased, and a signif- 
icant sub-gap current is visible. The difference between 
the two characteristics, which represents the contribu- 
tion of the atomic contact to the dissipative current, is 
shown in the inset with green open symbols. It is fit- 
ted with the theory 15-17 of Multiple Andreev Reflections 
(MAR) 18 in order to obtain the transmission of its con- 
duction channels 19 . The region \V\ < 200 /ieV, where 
sub-gap current was already found in the I — V charac- 
teristic of the Josephson junction by itself, was excluded 
from the fit. As a consequence, the accuracy on the de- 
termination of the transmissions of the channels is not as 
good as in experiments with atomic contacts alone 20 . It 
is here of the order of 1% for the largest transmission, and 
3% for the second largest. As an example, the contact 
corresponding to Fig. 5, called AC0 in the following, was 
found to have 3 channels, as most one-atom Aluminum 



contacts do; their transmissions, deduced from the fit of 
the I — V characteristic, are {0.95,0.45,0.10}. The cor- 
responding fit is shown in the inset of Fig. 5 as a solid 
line. 



IV. SWITCHING CURRENT MEASUREMENTS 

The supercurrent associated with the phase- 
dependence of the Andreev bound states energy is 
accessed through measurements of the "switching cur- 
rent" of the SQUID (see Fig. 5). The geometrical 
inductance of the loop, of the order of 20 pH 11 , can be 
neglected when compared to the Josephson inductances 
of the junction (Ljj = (fa/Io ~ 0.6 nH) and of the 
atomic contact (Lac ~ ^q/tA ~ 10 nH). Therefore, 
the phase 7 across the Josephson junction, the phase 8 
across the atomic contact and the flux phase ip — 2ir(f>/(f)o 
are related through ip = 8 — 7. 

The modulation of the switching current of the SQUID 
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In the experiments presented here, switching for the junc- 
tion alone occurs around 0.88, and 7 (s^) « 0.37T. 
As a consequence, switching for the SQUID occurs when 
I - Iat.c. (<p + 7( s ™)) ~ I s'd : tne average switching cur- 
rent is that of the junction alone, and the modulation cor- 
responds, in a first approximation, to the current-phase 
relation of the atomic contact, with a phase offset. The 
modulation of the switching current by the flux is there- 
fore a direct measurement of the current-phase relation of 
the atomic contact and not, as one could believe at first 
sight, a measure of the switching current of the atomic 
contact alone. 

Initially we performed switching measurement using a 
standard pulse technique: a train of current pulses of 
normalized height s — I /la is applied through the bias 
line, while monitoring the voltage across the SQUID to 
reveal switching. The total number of voltage pulses for 
a given current pulse train is recorded with a counter, see 
Fig. 6. In order to have detectable voltage pulses even if 
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Figure 5: Top: Current-voltage (J — V) characteristics of the 
Josephson junction (black open symbols) and of a SQUID 
with an atomic contact labeled AC0 (red filled symbols). The 
switching current of the SQUID, indicated with an arrow, 
is modulated by the applied flux. Bottom: same data, and 
(green open symbols) I — V characteristic of the atomic con- 
tact alone obtained by subtraction of the one of the junction 
from that of the SQUID. The transmissions {0.95,0.45,0.10} 
are found by fitting this I — V characteristic with the theory 
of Multiple Andreev Reflexions (green solid curve) , excluding 
the region \V\ < 200 ueV where the I — V characteristic of the 
Josephson junction presents resonances. 



can be qualitatively understood as follows 11 : when a flux 
is applied, the flux phase ip drops mainly across the weak- 
est link in the loop, i.e. the atomic contact. In contrast, 
when a bias current / is applied, it flows essentially in 
the arm with the largest critical current, which is the 
one with the Josephson junction. As a consequence, 7 ~ 
7(s) = arcsin (s) and 6 ~ ip + j(s), with s = I/Iq the nor- 
malized current and Iq the critical current of the Joseph- 
son junction alone. The current circulating in the SQUID 
loop, which is due to the phase-biased atomic contact, 
is therefore I at . c . (<p + j(s)) = ^ (kI A (ip + 7(5), n), 
and the total current through the Josephson junction 
I — Iat.c. (</? + 7(s)). Since Ia <C Iq, the switching cur- 
rent of the SQUID will be close to the switching current 
I'sw = s id^o of the junction alone, and 7(5) ~ 7 (s^) . 



Figure 6: Standard switching probability measurement 
scheme: a train of N pulses with normalized height sIq and 
duration t p are applied on the bias of the SQUID. The repeti- 
tion rate is 1/T r . The 40% lower plateau following each pulse 
holds the voltage at a finite value if switching has occurred, 
hence facilitating detection. The number n of voltage pulses 
resulting from switching events is recorded by a counter de- 
tecting crossings through a threshold value (dotted line). The 
switching probability is then P sw = ft . 

switching occurs at the end of the bias pulse, the actual 
bias pulse is followed by a 40% lower plateau that "holds" 
the voltage at a finite value for a time long enough for it 
to be detected. The height of this plateau is such that 
the switching probability during its duration is negligi- 
ble. We used a bias pulse duration t p =1 us, and 5 us- long 
"hold" plateaus. In general, the repetition period T r was 
20 us. The switching probability P sw (s) is obtained as 
the ratio of the number of measured voltage pulses to the 
number of bias pulses (typically 10 4 ). The inset of Fig. 7 
shows the resulting switching probability P s func- 
tion of the normalized pulse height s, for two values of 
the flux phase corresponding to switching currents close 
to its maximum and minimum values: (fx = 0.207T (black 
symbols) and if2 — 0.977T (orange symbols). Whereas 
the shape of the former (at (pi) corresponds to what is 
usually found on Josephson junctions, with a sharp vari- 
ation of the probability from to 1 as the pulse height 
is increased, the latter (at ^2) presents an unusual long 
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Figure 7: Switching probability P aw measured using the stan- 
dard method (see Fig. 6) as a function of flux phase tp and 
of the normalized pulse height s for the SQUID with contact 
ACO ({0.95,0.45,0.10}). Inset: P sw (s) for fluxes tpi = 0.20tt 
(black, open symbols) and tp <2 = 0.97-7T (orange, full symbols). 
Solid curves correspond to theoretical predictions. 



foot. The two types of behaviors are seen in the main 
panel of Fig. 7, a color plot of P sw with tp and s on the 
axes. Whereas the transition from P sw ~ (purple) to 
P sw ~ 1 (red) is abrupt for most values of tp, an anoma- 
lous region with an almost constant intermediate step at 
P sw ~ 0.1 (blue) is observed for 0.7ir < tp < I.Itt. In 
this region, the analysis of the histogram of time delays 
between switching events shows that they are correlated, 
with bunches of switching events and long "blind" periods 
without any switching. This behavior is reminiscent of 
the "blinking" observed in the fluorescence of molecules 
or quantum dots 21 ' 22 . 

Correlations disappear when a strong bias pulse that 
forces the system to switch (pulse height 30% higher that 
the measurement pulse, see Fig. 8) is applied before each 
measurement pulse (of course, the corresponding forced 
switching events are ignored in the counting). During 
the prepulses, a voltage develops across the atomic con- 
tact, and both the energy and the occupation of the An- 
dreev levels evolve at the Josephson frequency. When 
the current is reset to zero, the Andreev levels are back 
to their energies ±£U (S, Ti) , but the system is left in an 
out-of-equilibrium situation, which on average is always 
the same. The measurement pulse is applied a time At 
after this reset. Using such "switching prepulses", the 
switching probability is independent of T r , and the his- 
tograms of time intervals between switching events do 
correspond to independent processes. As reported in the 
main text, the corresponding curves P sw (s) then display 
a well defined intermediate plateau, as seen for flux tpi 
in the inset of Fig. 9. 
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Figure 8: Measurement scheme used in the following and for 
the data shown in the Letter. Short prepulses 30% higher 
than the actual measurement pulses cause the system to al- 
ways switch. With such pulses, successive switching events 
are uncorrelated. 



V. CALCULATION OF SWITCHING 
PROBABILITY 

We now give the details of the theoretical model used 
in the main paper to describe the switching probability 
P sw (s, tp) - It is related to the switching rate T(s, tp) and to 
the pulse duration t p by P sw = 1 — exp(— Tt p ). The phase 
7 across the Josephson junction is a dynamical variable 
governed by a Langevin equation, equivalent to the one 
obeyed by the position of a massive particle evolving in 
a "tilted washboard potential" 23-25 . The total potential 
of the SQUID is given by 9 : 



U(j) = -Ej cos 7 - Ejsj + ^2 c i E A{l - 



(6) 



where the first term is the Josephson energy of the tunnel 
junction, with Ej = I^tpo, the second one is the energy 
arising from the coupling to the bias source, and the last 
term is the total Josephson coupling introduced by the 
atomic contact, which depends on the configuration of 
the Andreev levels (cj = — 1 if channel i is in its ground 
state, c, = 1 in excited singlet, and Cj = in an odd con- 
figuration). Since Ej/A ~ 5.7 ^> 1, the two first terms in 
Eq. (6) dominate, and the shape of the tilted potential re- 
sembles that of a single Josephson junction, with slightly 
modified barrier height AU(s, tp) and plasma frequency 
u) p (s, tp). As in Ref. 9, switching can be fitted with 
thermal escape theory 26 , with a rate T = A exp(-B) 
with B = AU /ksT and A ~ oj p /2tt. When the atomic 
contact is open, the plasma frequency in the tilted po- 
tential is lu p — ojq(1 — s 2 ) 1 / 4 , with uq ss yloJ^poCjj, 
Cjj = 0.21 pF the Josephson junction capacitance esti- 
mated from its area 11 , and AU ~ (4y/2/S)Ej(l ~ s) 3/2 
the barrier height 27 . Fitting P sw (s) for the Josephson 
junction alone gives Io — 553. 7nA (luo/2tt 15 GHz), 
and an effective temperature T = 100 mK 28 . When an 
atomic contact is formed, precise comparison between 
experiment and theory is performed using a numerical 
determination of the barrier height, and of the plasma 
frequency ui p (s, tp) from the semi-classical calculation of 
the energy levels in the actual potential. More details on 
the SQUID potential are given in Appendix A. 

The predictions for the P 3W (s) curves of the SQUID 
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Figure 9: Color plot of the switching probability P 3W (s,tp) 
using switching prepulses (see Fig. 8, At = 0.5 us) for the 
SQUID with contact ACO with transmissions {0.95,0.45,0.10}. 
The color scale is the same as for Fig. 7. The white curves 
show the lines corresponding theoretically to P BW — 0.5 for 
the pristine contact ({0.95,0.45,0.10}, solid line) and for the 
contact with the first channel poisoned (dashed line). In- 
set: cut at tp = 0.97-7T (black line in the main panel). 
Black symbols: P S w{s) as measured; solid curves: theoreti- 
cal curves, Ps W {s) (leftmost, blue online) for the pristine con- 
tact ({0.95,0.45,0.10}) and Ps W (s) (rightmost, green online) 
for the poisoned contact, i.e. without the contribution of the 
most transmitting channel; intermediate line (black): fit of 
the data with the linear combination given by Eq. (7) with 
p = 0.36. 



with contact ACO are shown as solid lines in the insets 
of Figs. 7 and 9. Whereas the curve taken at tp = tp 1 = 
0.207T in Fig. 7 is well fitted by theory, the curve taken 
at tp = if2 = 0.977T has in common with theory only 
the value of s where the switching probability starts to 
rise. A central point of the paper is that, when pre- 
pulses are used (Figs. 9), the curve with a plateau can 
be very precisely accounted for by the weighted sum of 
two P sw (s) curves corresponding to two different config- 
urations of the contact, as shown in the inset of Fig. 9: 
P sw (s) = (1 - p)P° w (s) + pPL(s) , with p = 0.36, as 
well as in the inset of Fig. 2 in the Letter. The first one, 
P® w (s), is the one predicted for the pristine contact; the 
second one, P s 1 u; (s), is that of the poisoned contact, i.e. 
with its most transmitting channel in an odd configura- 
tion. This is the case in the whole flux region where the 
switching curves have an intermediate plateau, as shown 
in the main panel of Fig. 9. As in the Fig. 2 of the Letter, 
we compare the data with the two lines corresponding 
to the equations P s °^ 1 (s, tp) = 0.5. In the regions where 
P sw (s) has a standard shape, P sw (s)=0.5 occurs at the 



position predicted for the pristine contact (P° w (s)). In 
the region 0.77T < tp < 1.1%, P sw (s) has two steps (which 
appear in the figure as color gradients): one occurs at 
the position where P^ w (s, tp) = 0.5 (solid line), the other 
one at the position where P^ w (s,<p) = 0.5 (dashed line), 
showing that 

P sw (s, <p) = (1 - P &))PL(s, v)+P (<P) PL («, V>) • (?) 

In Fig. 10, similar comparisons are presented for three 
other SQUIDs formed with atomic contacts having one 
channel almost perfectly transmitting and the other ones 
with transmissions lower than 0.7. There again, in a 
broad phase range around 7T, P sw (s) shows a plateau de- 
limited by the predictions for the pristine (solid line) and 
for the poisoned configurations, i.e. with the more trans- 
mitting contact in an odd configuration (dashed line) . 



VI. MEASUREMENTS OF POISONING 
DYNAMICS 

A. Method 

The previous experiments demonstrate that, in a cer- 
tain parameter region, the system has a finite probability 
p to trap a quasiparticle and end in an odd configura- 
tion after the prepulse. In the following, we describe 
the experiments exploring the dynamics of trapping and 
untrapping, by varying the delay At between the pre- 
pulse and the measurement pulse. Data illustrative for 
the method, taken on a contact with transmissions {0.91, 
0.62, 0.15}, are shown Fig. 11: in a flux region exhibiting 
poisoning, the left panel shows P sw (s) for a short and a 
long delay At, and the right one the complete dependence 
of P sw on At at a bias value s = 0.826, for which P® w w 1 
and Pg W ~ 0. The data are well fitted with an exponen- 
tial dependence P sw ( At) = P oo +(P -P oo )exp{—At/T 1 ). 
One can then extract the initial poisoning just after the 
prepulse p = 1 — P Q , the asymptotic value at long times 
Poo = 1 — Pqo , and the relaxation time T\ . 

As discussed in the Letter, the dynamics of quasipar- 
ticle poisoning is strongly phase dependent. Measur- 
ing this phase dependence is however not trivial, be- 
cause during the pulse sequence shown in Fig. 8, 5 is 
not constant: it takes a value close to tp during At, then 
reaches ~ tp + 7 ( s iw) during the measurement pulse 
(7 { s sw) — 0.37r). As shown below, relaxation is very 
fast when S > I.Stt, so that data taken at tp > tt are 
dominated by the relaxation at measurement. 

To correctly measure the phase dependence of the re- 
laxation, we have therefore elaborated a different proto- 
col, shown in Fig. 12. We set the flux imposed by the 
external coil to a value tp i such that po is far from 0, and 
that the relaxation times at tp^ and at tp i + 7 (s^J are 
long, so that relaxation during the measurement can be 
neglected. We first take a switching curve P sw (s) and 
fit it with a weighted sum of two shifted curves as in 
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Figure 10: Color plot of the switching probability P sw (s,ip) using switching prepulses (see Fig. 8, with At = 2 /is) for SQUIDs 
made with three different atomic contacts, each having one channel with a transmission close to 1: AC1, {0.994,0.10,0.10} (same 
data as in the Letter); AC2: {0.998,0.59,0.13} and AC3: {0.993,0.65,0.46,0.10,0.05}. The white curves correspond theoretically 
to P sw = 0.5 for the pristine contact (solid line) and for the poisoned contact (dashed line). 
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Figure 11: Left panel: Switching probability of contact {0.91, 
0.62, 0.15} as a function of bias pulse height at fixed ip, for 
short (At = lps, green solid squares) and long (At = 580 us, 
red open circles) delay between the prepulse and the mea- 
surement pulse. Right panel: Dots: Evolution with At of 
the plateau height, measured at s — 0.826 (vertical line on 
the left panel). Solid line: exponential fit with P sw (At) = 
Pov + (Po - Poo) exp (— At/Ti) . 

the inset of Fig. 9. We then fix a working point sq cor- 
responding to the intermediate plateau, and determine 
the switching probabilities corresponding to the pristine 
contact P® w (sq) and to the poisoned contact P~ w (sq) . 
For all the data presented in the following, the work- 
ing point was chosen such that Pf? w (sq) = 1.00 and 
0.01 < P} w (s ) < 0.12 (except for the data taken at 
higher temperatures, where the rounding of the curves 
becomes comparable to their width, and P^ w (sq) can be 
as high as 0.32). The probability p to be in an odd con- 
figuration is then inferred from 

PsM = (i - p)p°M +pPl (so) , (8) 

which is slighlty more precise than the identification of 
p with 1 — P sw {so)- A measurement similar to that pre- 
sented in Fig. 11 is then used to characterize the system 



at flux ipi : 

p( i p h At) = £^?(p (tp i ),At) (9) 

where we have introduced the function (po , t) ac- 
counting for an exponential variation starting from po 
and during time t, with characteristic time T\ and asymp- 
totic value poo : 

(po, t)= Poo + (po - Poo) exp {-t/Tx) . (10) 

Using a dc flux pulse applied through the fast flux line 
between the prepulse and the measurement pulse, the 
flux phase is changed to a value ip w for a time t w (in 
practice, we leave a 1 us delay between the prepulse and 
the flux pulse to let the system stabilize; we also leave 
a delay of t° =15 us between the flux pulse and the 
measurement pulse to get rid of ringing effects after the 
fast flux pulse). Since the bias current is during this 
dc flux pulse, the phase across the atomic contact is then 
simply 5 = (f w 29 . 

From the switching probability P sw (t w + t° ) we cal- 
culate using Eq. (8) the probability p (t w + t°,) to be in 
an odd configuration after the complete sequence. We 
observe that p (t w + 1^) is not an exponential function 
of t w . The reason is that this probability results from 
two exponential relaxations with different parameters, as 
illustrated at the bottom of Fig. 12. The initial value 
p(0) = po {'■Pi) results from the prepulse applied at phase 
flux tfi (the function po (tp^) is discussed farther). Fol- 
lows an exponential evolution at phase flux tp w during 
t w , leading to 

p(t w )=S^ ) ) ( P (0),t w ). (11) 
During the last i° = 15 us, the phase flux is ipi, so that 
P {t w + t° w ) =^°&)(P ( 12 ) 
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Figure 12: Pulse sequence used to measure the phase depen- 
dence of the relaxation process. The first line corresponds 
to the signal applied to the current bias, the second one to 
the fast flux line. The two last lines sketch the correspond- 
ing evolution of the phase 5 across the atomic contact and 
of the poisoning probability p(t). The prepulse causes sys- 
tematic switching, and the "running away" of the phase is not 
represented. The dotted line in S(t) during the measurement 
pulse indicates that either the SQUID switches and the phase 
runs away, or it does not switch and S simply follows the bias 
current. The prepulse and the measurement pulse are always 
applied at the same flux phase tfii. The probability p (t) starts 
from p (0) after the prepulse, a value that depends on the flux 
<Pi; a flux phase ip w is then applied during a time t w , and p(t) 
relaxes exponentially with a time constant Ti (ip w ) towards 
Poo {<Pw), reaching p (t w ) . In the last step, the phase flux is re- 
stored to ifi during t%, = 15 us and p (t) evolves with the time 
constant Ti (ipi) towards poo (y>i), reaching finally p(t w 
the actual value accessed by the measurement (the schematic 
corresponds to a situation where p (t w ) < p^ (ifi), hence p (t) 
increases in the last step). In the time interval between the 
prepulse and the measurement pulse, the bias current h is 
zero, and the phase 5 across the atomic contact is equal to 
the flux phase: 5 = ip w during t w , then 8 = tpi during t{j,. 



Since the parameters p(0), T\ (cpj) and p^ (ifi) have 
been determined in the first measurement without the 
flux pulse, Eq. (12) can be used to deduce p(t w ) from 
p (t w + t° ) . The function p(t w ) is then an exponential, 
as expected, and its fit with Eq. (11) yields the asymp- 
totic poisoning p^ (6) and the relaxation time T\ (<5) . 

Using this more complex procedure that involves the 
fast flux line to apply dc flux pulses, the relaxation can 
be probed with the prepulse and the measuring pulse 
both acting always at the same flux. The initial value 
of the poisoning probability is therefore always the same, 
and the working point can be chosen such that relaxation 
during the measurement pulse plays no role (note that, as 
long as this last criterium is observed, p^ (<5) and T\ (8) 
do not depend on the working point). 



B. Temperature dependence 

The functions poo (5) and Ti {$) at various tempera- 
tures, are shown in Fig. 13 for contact AC1. Whereas 
Paa {$) hardly changes, T± (8) drops rapidely with tem- 
perature, and relaxation is almost instantaneous above 
220 mK. 
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Figure 13: Relaxation time Ti and asymptotic poisoning 
probability p^, as a function of the phase S — <p w applied 
during the dc flux pulse imposed for a time t w between the 
prepulse and the measurement pulse. The data are taken on 
the atomic contact AC1 ({0.994,0.10,0.10}), and at four dif- 
ferent temperatures T indicated in the figures. 



C. Revisiting P sw (s,ip) data 

Coming back to the results presented in the left panel 
of Fig. 10, the boundaries of the phase range in which poi- 
soning occurs can be understood in the light of Fig. 13: 
the left boundary corresponds to the phase cp^ at which 
7i (ipl) becomes comparable with At = 2 /xs, leading to 
significant relaxation in the time interval between the 
prepulse and the measurement pulse. The nature of 
the right boundary is different: it corresponds to the 
phase ipR at which Ti (7 (s) + <pn) becomes comparable 
with t p = 1 fis (we recall that 7 (s) is the phase across 
the Josephson junction during the measurement pulse), 
leading to significant relaxation during the measurement 
pulse. In contrast with the left boundary, the position of 
the right one depends slightly on s through 7 (s) , which 
explains that it is slightly tilted. Hence, in the simplest 
procedure where the flux is the same during the whole se- 
quence, the effect of relaxation during the measurement 
pulse becomes predominent for ip > it, explaining why 
the intervals on which poinsoning is observed in Fig. 10 
is not centered at ir. 
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Figure 14: Relaxation time T\ and asymptotic poisoning 
probability p^, measured at T = 20 mK, for five different 
atomic contacts with the shown transmissions. 



D. Revisiting data without prepulses 

The anomalous statistics of the time intervals between 
switching events mentioned above finds also an expla- 
nation. The data presented in Fig. 7 are taken without 
prepulse, so that if poisoning occurs after a measurement 
pulse, the system remains poisoned during a time T\ in 
average, and switching is suppressed, giving rise to long 
"blind" periods. When the system unpoisons, switch- 
ing occurs again, and, the probability to get poisoned 
again being rather small, several switching events occur 
in a row. When the repetition period T r is increased, 
the probability to escape from a poisoned configuration 
before the next measurement increases. And since the 
switching rate is larger when the system is not poisoned, 
the average switching probability increases, as observed. 



E. Raw data on several contacts 

Similar data taken on a variety of atomic contacts are 
shown in Fig. 14. The phase interval in which poisoning 
occurs reduces when the transmission of the most trans- 
mitting channel diminishes. The same data are replotted 
as a function of in Fig. 4 of the Letter. 



F. Multiple poisoning 

In a contact with more than one well transmitting 
channel, poisoning can affect several channels at once, as 
shown in Fig. 15 where the switching probability presents 
two intermediate plateaus. The first one near 0.5 corre- 
sponding to either one of the two first channels being 
poisoned (they have very similar transmissions), the sec- 
ond one at 0.95 to the situation in which both are. 
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Figure 15: Measured switching probability as a function of s 
for a contact with two well transmitting channels: {0.96, 0.95, 
0.60, 0.34, 0.30, 0.29, 0.27, 0.26, 0.24, 0.2}, taken at a flux 
within the poisoned region. The first plateau (P sw — 0.55) 
can be attributed to situations where one of the two channels 
with transmission ~ 0.95 is poisoned, while the second one 
corresponds to situations where both are poisoned. 



G. Initial poisoning 

We have measured the initial poisoning po as a func- 
tion of the phase 8 across the atomic contact AC1 
({0.994,0.10,0.10}). The measurement protocol is shown 
in Fig. 16: the flux phase ip pp applied till the end of the 
prepulse is varied, then the flux is reset to ipi. Here again, 
the effect of relaxation at tp^ during the 15 us between 
the end of the flux pulse and the measurement pulse has 
been corrected for. The result of the measurement is dis- 
played in Fig. 17 (for the chosen value of ifi, indicated 
with a dashed line, the relaxation before measurement 
is characterized by T\ — 167 us and Poo = 0.12). This 
irregular pattern is responsible for the vertical stripes in 
the data shown in Fig. 10. Similar patterns, with slight 
dependence on the duration of the prepulse and on the 
rate of decay of the current at the end of the prepulse, 
were found on the other contacts, but the influence of the 
different parameters could not be deconvolved from the 
data. 



VII. QUASIPARTICLE POISONING ANTIDOTE 

The experimental data indicate that when the phase 
6 across the atomic contact is driven far from tt, the 
system quickly relaxes to the ground state. Based on 
this result, we have developed a procedure to get rid 
of trapped quasiparticles. The principle, illustrated in 
Fig. 18, is simply to sweep the flux phase over 2ir, with 
a pair of symmetric "antidote" pulses with amplitude tt 
each, hence insuring that whatever the starting phase, 
d crosses or 2-7T. Poisoning is cured by the antidote 
pulses, as shown in Fig. 19 where we compare relaxation 
curves with and without them at ip w = n for contact AC1 
({0.994,0.10,0.10}), the flux phase value where the poi- 
soning probability was the highest. Whereas data taken 
as in Fig. 11 (right panel) show initial poisoning with a 
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Figure 16: Pulse sequence used to measure the phase depen- 
dence of the initial poisoning probability po- The first line 
corresponds to the bias current, and the second one to the 
fast flux line. The third line sketches the evolution of the 
phase 8 across the atomic contact, and the last one that of 
the poisoning probability. The quantity that is studied is the 
poisoning probability po (<p PP ) just after the prepulse, which 
depends on the phase flux ip pp applied to the contact during 
the prepulse. Note that the measuring pulse gives access to 
p (t^) , which includes relaxation at flux ipi during the last 
t^, — 15 us. This relaxation is taken into account to obtain 
Po (<Ppp)- 
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Figure 17: Initial poisoning po as a function of the phase 
8 = ifi pp across the atomic contact AC2 at the end of the 
prepulse. The pulse sequence used to gather this data-set is 
shown in Fig. 16. The dashed line indicates the value of the 
phase S across the atomic contact when only the dc flux ipi is 
applied. 
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Figure 18: Pulse sequence used to measure the phase de- 
pendence of the relaxation process after the application of 
"antidote pulses". The first line corresponds to the current 
bias, and the second one to the fast flux line. The third line 
sketches the evolution of the phase 8 across the atomic con- 
tact, and the last one that of the poisoning probability p (t). 
The antidote pulses are pairs of dc flux pulses applied in both 
directions, with an amplitude corresponding to a phase ex- 
cursion of 7r on the phase 8. Their total duration was 0.5 us; 
they are applied 1.5 us after the prepulse. During the pair 
of antidote pulses, the phase visits regions where a trapped 
quasiparticle escapes rapidly, and the poisoning probability 
decays to 0. 



ity was checked on many atomic contacts. These data 
prove the robustness of the exponential behavior, and 
that poisoning associated with residual quasiparticles has 
the same effect whatever the initial configuration. 

Antidote pulses were employed to record P sw (s,<p) on 
the SQUID with contact AC0, as shown in Fig. 20. We 
used the simplest current bias pulses, shown in Fig. 6, 
with antidote pulses applied on the flux line 1 us before 
the measurement pulse. The region displaying poisoning 
in Fig. 7 has completely disappeared, and switching is 
regular for all values of the flux phase ip. 

The procedure presented here allows to obtain with 
very high probability unpoisoned configurations, at least 
just after the antidote pulses. Nevertheless, attempts to 
induce transitions from state |0) to state |2) with mi- 
crowaves applied on the fast flux line were unsuccessful 
in this experiment. 



probability po — 0.37, the poisoning probability extrap- 
olates to at t w — when antidote pulses are applied 
(Fig. 19, top curve), indicating that poisoning is absent 
just after the pulses. Remarkably, the subsequent evolu- 
tion is identical in both situations, tending exponentially 
to p = 0.23 with a time constant T\ — 220 us. The iden- 
tical time constants and asymptotic poisoning probabil- 



VIII. MECHANISM FOR INITIAL POISONING 

The system always switches during the 0.1 us- long pre- 
pulse, and a voltage of the order of V = 0.1 mV devel- 
ops across the atomic contact. The phase S then varies 
rapidly, at the pace V/tpo, and multiple Andreev reflec- 
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Figure 19: Symbols: Poisoning probability as a function of 
time t w spent at phase n, for contact AC1 ({1,0.07,0.07}). 
Top curve is taken with the pulse sequence shown in Fig. 12, 
and shows frequent poisoning (p (0) = 0.37). Bottom curve 
has "antidote" flux pulses after the current prepulse (see 
Fig. 18), which suppress initial poisoning (p (0) = 0). Solid 
lines are exponential curves with identical time constants 
Ti = 220 us and asymptotic value p = 0.23. 
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Figure 20: Color plot of the switching probability P 3W (s,ip) 
using an antidote flux pulse 1 us before the measurement 
pulse, and without prepulse, for the SQUID with contact AC0 
with transmissions {0.95,0.45,0.10}. The color scale is the 
same as for Fig. 7. The white curve show the lines corre- 
sponding theoretically to P 3m — 0.5 for the full contact. This 
plot is to be compared with Fig. 7. 



tions carry a current across the atomic contact, which 
creates quasiparticles on both sides of the contact 30 . For 
a contact of unit transmission, two quasiparticles are be- 
ing created for each turn of the phase. The resulting 
number of quasiparticles is therefore 2 x 0.1 us x 0.1 mV x 
2e/h ~ 10 4 . When the bias current is reset to zero, the 
phase stops in a local minimum of the global potential, 
and the configuration of the Andreev levels freezes. The 
experiment indicates that after this step, the population 



of the Andreev states generally does not correspond to 
the equilibrium situation: there is a poisoning probabil- 
ity po from which relaxation is subsequently observed. 
The exact process of poisoning after switching is not un- 
derstood presently. One can however qualitatively grasp 
the mechanisms that can be relevant using the picture of 
the "Andreev elevator" that describes transport in the fi- 
nite voltage state, at least for small voltages 30 : the phase 
varies linearly in time and the Andreev levels, in the semi- 
conductor representation, coincide periodically with the 
gap edges of the continuum (when S = 2nn) , loading the 
lower level with a quasiparticle at the lower continuum, 
and unloading it in the upper one. The quasiparticle 
is sometimes transferred from the lower level to the up- 
per one by Landau-Zener tunneling. In this scenario, 
the system is always in an even configuration. However, 
many quasiparticles are created in the continuum, and 
the "loading" (or "unloading") process can fail, leading 
to a situation where, after passing 6 — 2wr, the lowest 
state is unoccupied while the upper one has emptied (or 
vice versa), leading to an odd configuration. If the sys- 
tem returns to the zero voltage state at this stage, the 
system stays trapped in the odd configuration. The os- 
cillating structure for po ((p pp ) in Fig. 17 suggests that 
interferences during the phase dynamics may also play a 
role. Further work is clearly needed to clarify this phe- 
nomenon. 



IX. MODEL FOR POISONING DYNAMICS 

We now discuss with the help of Fig. 21, the details 
of the model presented in Fig. 4 of the Letter. Quasipar- 
ticles can jump between the Andreev states and quasi- 
particle states in the continuum. The processes are in 
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Figure 21: Model for the dynamics of the population of the 
Andreev configurations. Left: in the excitation representa- 
tion. Right: in the configurations space, where the probabil- 
ity of the 4 configurations are given in itallic. Quasiparticles 
can jump in and out the Andreev levels at rates Tin and r out . 
If both levels of the Andreev doublet are occupied, the system 
can decay directly to the ground state (rate T2o). 

principle rather slow because they either require energy 
absorption (rate r out ), or the presence of quasiparticles 
in the continuum (rate Ti n ). Because the state with a 
double excitation is never observed, we assume that the 
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relaxation rate r 20 to the ground state is much larger 
than Tin and r out . This hypothesis is discussed in Ap- 
pendix B. 

The master equation for the populations of the ground 
state |0), 1 — p, and of the odd configurations |1 1) and 
|1 4-)i f/2 each, is (see right hand side of Fig. 21): 

d (1 d ~ P) =-2(1- p) T ia + p (Tout + r in ), (13) 

the term +pT ln corresponding to transitions from the odd 
configurations through the excited state |2) which relaxes 
very fast to the ground state. Assuming that Ti n and r out 
arc time- independent, which is necessary to obtain the 
observed exponential dependence of p(t), one gets equa- 
tion (2) of the Letter. The assumption that r; n , which 
is proportional to the quasiparticle density, is constant, 
means that quasiparticles created in the voltage state af- 
ter switching have disappeared at the time scales probed 
by relaxation experiments. Experiments with a succes- 
sion of several prepulses indicate that this is the case. 
One possibility is that they have already recombined. In 
the next paragraph, we argue that they have probably 
already diffuse away. 

To evaluate the decay of the density n qp of quasi- 
particles occurring due to just diffusion, we start from 
the number of quasiparticles created during the pre- 
pulse, 10 4 , as evaluated at the beginning of Section VIII. 
Within the duration of the prepulse, they spread over 
~ \]D qv x 0.1 us > ^D N x 0.1 us ~ 50 pm (with D qp = 

Dn/ (l — (A/E) 2 ^j the diffusion constant for quasiparti- 
cles of energy E in the superconducting state 31 , Dm the 
value in the normal state), i.e. over the whole area of 
SQUID (see Fig. 3(c)). The density at the end of the 
prepulse is therefore n qp rts 100 ^m~ 3 . Diffusion is then 
one-dimensional because it is limited by the 3 thin con- 
nections (inductive lines), and a rough estimate (using 
again the diffusion constant Dn instead of D qp ) indi- 
cates that the density decays within a few microseconds 
below 10 /im~ 3 , which is the typical background value 
found for Al films at low temperatures in Refs. [32,33]. 
The quasiparticles created by the prepulse can therefore 
be neglected when describing the dynamics of poisoning 
relaxation. 



X. QUASIPARTICLE TRAPS 



that the diffusion time of a quasiparticle through it ex- 
ceeds the relaxation time of its energy. 

On the other hand, the decoupling from dissipative 
parts of the circuit is an important requirement 34 in or- 
der to explore the physics of the Andreev qubit 35 ~ 37 , and 
the use of normal metal traps is probably not advisable. 
As a consequence, there would always be some quasipar- 
ticles ready to get trapped in the Andreev bound states. 
One method to get rid of them is to use the "antidote 
pulses" described in Section VII. Another method is to 
have connecting wires with a lower gap than the SQUID. 
We checked experimentally, using bilayers of Al and Cu, 
that it allows to get rid of poisoning. 



Appendix A: SQUID potentials 

The SQUID potential is given by Eq. (6). If transi- 
tions between configurations are slow compared to the 
plasma frequency, which is the characteristic timescale 
for the dynamics in U^), the escape rate is deter- 
mined by the potential corresponding to the instanta- 
neous configuration 38 . As an example, we plot in Fig. 22 
the three potentials for a contact with a single channel 
with transmission 0.99, at ip/n = 0.9, and for s = and 
s = 0.87. The plot at s — matters for the evaluation 
of the transmission rates during At, whereas the plot at 
s = 0.87 aims at illustrating why the switching rates are 
very different in the ground state of the channel and in 
the odd configurations. The energies Ej = 1.14 meV 
and A = 0.193 meV correspond to the parameters of the 
Josephson junction of the experiment. The distance in 
energy between the potentials is the smallest when S = ir, 
i.e. 7 = 7r — ip = O.ln. We also indicate the ground state 
energy in each potential, as obtained from semi-classical 
calculations. The quantum fluctutations of 7 (and hence 
of S) are of the order of O.I71-. When the flux phase tp is 
different from and tt, as in the example of Fig. 22, the 
three potentials do not have their minima at the same 
value of 7, an effect that is expected to modify the tran- 
sition rates between configurations 11 . 

At s = 0.87, the barrier heights are strongly reduced. 
The barrier is the smallest in the ground state, and 
switching occurs only if the system is in the ground state, 
and not in the odd configuration, as is the case in Fig. 10 
for ip/n — 0.9 between the dotted line and the solid line. 



Let us note that in previous experiments with atomic 
contacts in SQUIDs 9 ' 30 we didn't observe poisoning. In 
these works, the superconducting loop was directly con- 
nected to very large normal metal electrodes, which acted 
as efficient quasiparticle traps. In the present case, our 
SQUID is contacted through very narrow and thin su- 
perconducting lines, and clearly the small pieces of gold 
placed in direct contact with the present SQUID body 
didn't act as efficient quasiparticle traps. In fact, a nor- 
mal metal electrode can only act as a trap if it is so large 



Appendix B: relaxation between Andreev states 

The excited doublet |2) is never observed in our exper- 
iments. We have attributed this fact to a large value of 
the relaxation rate r 2 o to the ground state. Two contri- 
butions to this rate can be evaluated. The first one is the 
rate F^™ of the process in which a photon is emitted into 
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Figure 22: Potentials of a SQUID with a single channel 
atomic contact with transmission r = 0.99, at (p/ir = 0.9 
and s — (top panel) or s — 0.87 (bottom panel). The three 
curves correspond, from top to bottom, to the atomic contact 
in the ground state |0) (blue), in the odd configuration |1) 
(green), and in the excited state \2) (red). Horizontal lines 
indicate the energy of the first level |&j p in each potential, 
with the plasma frequency uj p /2tt calculated semi-classically 
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Note that this equation actually differs by the last factor 
from the one in Refs. [34], because we have used to de- 
scribe the Andreev system the hamiltonian proposed in 
Ref. [36] instead of the one introduced in Ref. [41] which 
does not enforces charge neutrality. For the actual circuit 
design (Fig. 3), the environment impedance Z t is: 



(u) = jCjjLU + 



1 



jLjjw jLuj + 



(15) 

where Ljj — (po/Io, Cjj are the Josephson junction 
inductance and capacitance, respectively. The eigen- 
modes of the two LC circuits, at 1/ (2wyLC) ~ 0.5 GHz 
and l/(2ny/LjjC,jj) ~ 14.5 GHz, couple to give two 
peaks near 0.4 GHz and 17 GHz; between 0.5 GHz and 
15.3 GHz and above 16.6 GHz, the impedance ReZ t is 



smaller than 1 Q. The resulting rate T^™, shown with 
solid lines in Fig. 23, exceeds 1 MHz near 5 — tt when 
the transmission is large enough. Note that these esti- 
mations depends crucially on the value of the losses in 
the capacitor, which are modelled by r = 0.5 il , a value 
determined from low- frequency (0.5 GHz) measurements. 
When the Andreev frequency 2_E/i//i is large, the losses 
could easily be very different, and the rate much larger. 

The second contribution to is the rate T^g of the 
process in which a phonon is emitted. An upper bound 
to this rate 40 has been calculated in Refs. [37,41] in the 
limit E A < A: 



20 



(1 



E A 



pit 
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(16) 



with r pf ^ (E) = K p h E 3 the bulk electron-phonon relax- 
ation rate at energy E. Experiments 42-44 on similar Al 
films give Kph — 3/iS _1 K^ 3 . The corresponding predic- 
tions are shown with dashed lines in Fig. 23 for trans- 
missions 0.80, 0.95 and 0.995. Only for r < 0.80 does 
the rate T^q exceed 1 MHz for phases around tt. There- 
fore, the large value of T2o in the phase interval where 
poisoning is observed (at the MHz scale, which is the up- 
per bound to the measured values of T in ^ out , see Fig. 4 
of Letter) can be attributed to losses in the electromag- 
netic environment of the atomic contact and, to a lesser 
amount, to phonon emission. 
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Figure 23: Calculated prediction for the relaxation rate 
from the excited singlet, for transmissions r = 0.8, 0.95 and 
0.995. Two contributions are plotted: rate for photon emis- 
sion into the electronic environment ri™ (solid lines) and up- 
per bound 40 for the phonon emission rate Tjq (dashed lines). 
For r = 0.995, three peaks are visible in F^o 1 , corresponding to 
the Andreev frequency 2Ea / h approaching the low frequency 
mode of the electromagnetic environment at 0.4 GHz (central, 
broad peak) or coinciding with the high frequency mode at 
17 GHz (lateral, sharp peaks). For smaller transmissions, only 
the effect of the high frequency mode is visible. 

Given the values calculated for T20, one could expect 
sharp spectroscopic signatures of the transition |0) — > |2) 
when microwaves are applied. However, no clear trace 
of this effect could be found in the experiment. This 
could be due to imperfections in the implementation of 
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the microwave environment of the SQUID, or to addi- tional relaxation channels. 
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